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Visual  Information 


When  we  see  something,  we  are  not  seeing 
an  object,  but  our  senses  capture  sensorial 
information  (luminance  contrasts),  which 
should  then  be  transformed  into  an 
intermediate-level  representation  of  what 
gives  the  perceived  entity  its  individuality. 


VISUAL  INFORMATION 


raw  visual  data 


a 

image  segmentation 

~ n 


low  -level  interpretation  of  picture  elements 


a 

object  recognition 

■wwwvlBwwww^ 

a 

high-level  image  interpretation 
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Low-level  information  is  typically 
about  the  spatial  relationships  among 
primitive,  two-dimensional  visual 
features  such  as  observed  shape,  texture, 
and  composition  variability  patterns. 

Intermediate  information  describes  the 
properties  arising  from  forms  of 
organization  of  the  low-level 
primitives,  and  may  include  descriptions 
of  the  three-dimensional  spatial 
relationship  (location)  among  visual 
properties. 

The  overall  explanatory  process  is  thus 
broken  down  into  the  extraction  of  a 
number  of  different  observable  physical 
properties  (low-level  analysis),  followed 
by  a final  decision  based  on  these 
properties  (high-level  analysis),  what 
implies  breaking  down  the  perception  of 
meaningful  visual  marks  into  different 
explanatory  stages 


Visual  Information 


Formally  speaking,  a surface  is  a boundary 
of  separation  between  two  phases. 

In  its  turn,  a phase  is  a homogenous  mass  of 
substance,  solid,  liquid  or  gas,  possessing  a 
well-defined  boundary. 

When  we  have  two  phases  in  mutual 
contact,  we  have  an  interface. 


Visual  Information 


The  surface  of  solids  plays  a significant  role 
to  discover  the  way  they  have  been 
produced  and  the  way  they  have  been  used. 

Surfaces  have  two  main  properties:  shape 
and  texture. 


What  is  Shape? 


Scientists  have  traditionally  assumed  that  there  is  a roughly  fixed  set  or  vocabulary  of 
“supposed”  descriptive  visual  regularities  shared  by  a single  population  of  objects,  which 
are  also  distinctive  enough. 

Scientists  believe  that  what  they  see  is  a “seed”,  a “bone”,  a “bowl”,  “a  knife”,  the  “wall  of 
a house”,  a “prince  burial”,  etc.,  and  they  can  distinguish  between  different  kinds  of 
“bowls”,  different  kinds  of  “prince  burials”,  and  so  on. 
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TABLE  PREHENSION 


□ prehension  appliquec  □ tiisclkmtiM  median 

□ prehension  tiree 


Type 

Localisation 

Assoc  ic  a 

Orientation 

- 

jn.se  en  bond  in 

- 

sous  la  levre 

- element  plaslique 

- verticals 

- 

anse  en  boudin  a are  cimtre 

- 

sur  la  levre 

element  prehension 

oblique 

* 

anse  en  boudin  a are  coude 

- 

a prise  direete  sur  la 

element  en  relief 

hnrizontalc 

- 

anse  en  boudin  a angle  droit 

levre 

decor  en  relief 

indeterminable 

- 

anse  en  ruban 

* 

sur  te  bond 

decor  encreux 

- 

anse  en  ruban  a arc  eintitf 

- 

sur  le  col 

indeterminable 

- 

anse  en  ruban  a arc  coude 

- 

sur  la  pause 

- 

anse  en  ruban  a angle  droit 

- 

en  haul  de  la  pause 

- 

anse  en  demi-bobine 

- 

au  milieu  de  la  pause 

- 

arisen  ad a$eia» 

- 

on  bos  do  la  panse 

- 

mamelon 

- 

a cheval,  Levre/bord 

E 

- 

mamelon  allonge 

- 

4 cheval  levre/col 

e 

bi 

- 

ma  melon  prism  at  ique 

- 

a cheval.  livre/panse 

- 

prehension  en  demi-bobine 

- 

a cheval  borti/panse 

- 

prehension  en  demL-cylindre 

- 

a cheval  col/panse 

- 

prehension  en  H 

- 

au-dessus  dc  la 

- 

prehension  en  X 

care  no 

prise  plate 

- 

sur  tj  caiene 

prise  plate  a devcloppcmcnl 

- 

au-dessous  de  la 

arrondi 

carene 

prise  plate  a developpement 

. 

mtiachee  au-dessus 

rectangulaire 

de  Sa  carene 

indeterminable 

- 

rattachee  au^dessous 

prehension  tubulaire  a 

de  la  carene 

developpement  arrondi 

- 

indeterminable 

prehension  tub-til  aire  a 
developpeme nt  rectangu  la  i re 

Perforation  affectant  la  prehension 

" notnbre  tie  perforation 

□ som-cutanee 


Orient  til  in  it  pcrfo/rcdpicnl 

Orientation  ptrfo/prthension 

E 

s 

R 

E£ 

verticale 

oblique 

horizontal? 

indeterminable 

transversale 

longitudinals 

indeterminable 

— T ,\|  D tj  S i f'1,1  F i_i  1“  H " L1  L Dt  i-  i i C U L i f-4  ,r. 
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Tipologie  Analytique 


Figure  Ui 

1 : retouche  simple,  piofcmdc,  lateral?  : racloir 
lateral. 

2 : retouche  simple,  marginal?,  lateral?  : racloir 
lateral, 

3 : m£me  type  de  retouche  mais  late  rale  et  trans^ 
versale  - racloir  latero-transversaL 

4 ; retouche  simple*  profondeL  transversal?  ; 
racloir  transversal. 

7 : meme  type  que  4,  mats  en  reiouche  margi- 
nalc  ; racloir  transversal. 


6 : retouche  simple,  profonde,  lai£rate.  Le  bard 
destre  «t  denticul£-  Un  front  £,  peine  cbauchc  ne 
peul  etre  interpret?  cornrne  grattolr  II  s’agk  d'un 
racloir  denticulc. 

5 : le  bord  dcxirc  est  un  racloir  lateral  pan  jet.  Le 
bord  seriestre  prfoeme  une  reiouche  abrupie,  mar- 
ginale,  qui  d Am  sa  panic  proximal?  rencontre  une 
cassure  concave  amenagfce  : il  s'agit  d'un  bee. 


Para  su  descripcion  tenemos  en  cuenta  los  siguientes  parametros: 
o Forma  de  la  pieza  (en  vision  frontal): 


o 


Triangular  V ; rectangular 
o Filo  (en  vision  frontal): 


I "■ 

; trapezoidal  V#-j 


Rectllineo 

Perfii  del  filo: 


>*>>■•  :-'.l 

\ • -’  >-  •7  ; 


curvo 


...  j 


S t met  r ico  f i 'f ; a s i met  ri  co  1 
o Base  (en  vision  frontal): 


Redondeada 

o Seccion: 


; recta  \L 


Rectangular 


; en  media  iuna  ^ 


; oval  ; circular 


12 


13 


BOrdS 


(rypologie  d'aprAft  Lefr*nc_  2007) 


This  way  of  identification-based 
explanation  seems  then  a tricky  way  of 
solving  any  archaeological  research 
problem. 

It  pretends  to  explain  what  has  been 
“seen”,  not  in  terms  of  their  visual 
characteristics,  but  in  terms  of 
subjective  recognition. 

Nevertheless,  what  we“see”  in  the  real 
world  are  not  stones,  walls,  pit  holes, 
mounds,  buildings,  pottery  sherds, 
plants,  animal  carcasses,  or  anything 
like  but  a hierarchized  organization  of 
visual  marks  and  higher  level  cues  to 

explanatory  categories.  127/ 


What  is  Shape? 


Shape  is  a quantitative  property 
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Quantitative  Properties 


QUANTITY:  any  property  that  can  be 
expressed  in  terms  of  an  intensity. 

QUALITY:  any  property  that  can  only  be 
expressed  in  terms  of  presence/absence 


Quantitative  Properties 


Verbal  language  is  not  appropriated  for 
describing  intensities: 

“a  lot  of”, 

“small” 

“ large 99 

“ tall ” 


Quantitative  Properties 


There  is  not  a “ mathematical ” reality , nor 
“ mathematical  concepts  *\ 

There  are  real  things,  empirical  properties 
that  SHOULD  be  described  in  terms  of 
their  intensities 


Measuring 


Measuring  is  the  operation  of  assigning 

numbers  that  represent  the  degree  or 
intensity  a property  is  present  in  some 

entity 


Measuring  Instruments 


It  is  a device,  a procedure  or  a set  of  logical 
formal  operations  that  compare  a reference 
ordination  to  the  actual  state  of  the  property 
at  the  measured  entity 


Measuring 


F ast  everything  observable  in  reality  is  a 

quantity. 

Happiness,  Love, 

they  are  always  refered  as  if  they  would  exist 

in“quantities” 

I am  “very”  happy. . . 

I love  you  less  and  less. . . .) 

But  we  do  not  know  how  to  measure  those  properties 


Measuring  Shape 


We  do  not  have  the  habitude  of  measuring 
shape  in  archaeology,  but  it  is  far  easier 
than  measuring  happiness  or  love 


Shape  is  not  Size 


One  accurate 
measurement  is  worth 

a thousand 
expert  opinions 

Grace  Hopper 


:Ssii 

onfine-behavtar 
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Ceramic  Head 

Stauffer-Allison 

collection 


Lyon  County 


Rad  ins 


Height 


; u < 

ii 

A partir  des  mensurations : 

ip,  ULf. prises  sur  les  d iff  ^rentes 

parties  dun  vase,  ont  calcutees 

les  27  variables : Pnses 

Lt  Lp  hp 

en  compte  pour  /‘analyse  factoriel/e. 


D 


2 H 


Fig.  6.3:  DiffSrents  exemples  de  mensurations 
d) : Ceramiques  (d'apres  Mohen,  1980)  c)  : Ceramiques  (d’apres  Whallon,  1982). 


Figure  1»  M^NumiU'iibi  Liken  on  Andra  Pradesh  vessels 
(mils  catfijtr 

enable  us  to  calculate  indexes  of  regularity  in  the 

i f • n .%  it  ii  .1  i ■ * 


Figure  2.  Measurements  taken  on  Spanish  vessels.  (1) 
Maximum  il  in  meter,  (2)  aperture,  (3)  netk  diameter,  (4) 
base  diameter,  155  total  height,  (rtf  maximum  diameter 
height,  (7)  lip  height,  (8)  lip  thickness,  (f  ) handle  width, 
(10)  handle  thickness,  (11)  height  of  the  decor  (distance 
between  the  bottom  and  ail  incised  groove),  and  (12) 
height  of  the  tone  which  has,  not  been  turned. 


Local  ion  oC 


L Maximal  Width  . I 

i1 


Fig.  I.  Conventional  positioning  and  basic  measurement  extraction  (from  Roe,  1964 
and  after  Sharon,  2007}. 


Figure  1.  Standard  fluted-point  morphometric  landmarks:  (mi)  maximum  length,  (mw)  maximum  width,  (Ig) 
extent  of  lateral  grinding,  (bw)  basal  width,  (bd)  basal  depth,  (fl)  flu  ting-scar  length,  and  (ml)  maximum  thickness. 


Variable  number 


8.  Maximum  length  (cm) 

<MAXLE> 

9.  Length  of  socket  (cm) 

<SOCLE> 

10.  Maximum  width(cm) 

<MAXWI> 

11.  Width  of  upper  socket  (cm) 

<UPSOC> 

12.  Width  of  lower  socket  (cm) 

<LOSOC> 

13.  Distance  between  maximum 
width  and  lower  socket  (cm) 

<MAWIT> 

14.  Weight  (g) 

< WEIGHT  > 


Nomenclatura  morfologica:  mKr_ 

a)  Apice  o punta.  b)  Aleron.  c)  Pedunculo  d)  Hombre- 
ra.  e)  Base  del  pedunculo.  f)  Cuello.  g)  Borde  lateral, 
h)  Borde  del  pedunculo,  (de  BAGOLINI) . 


Pin  ^nn 


Tipos  de  alerones,  segun  HUGOT 
a:  Agudos.  b:  Rectos,  c:  Obtusos 

Fig.  335 


Maximal  length 


Industria  osea.  Stordeur  1992 


< > 


thickness 

Di-sfeal  mvbcSlhi 


Maximal  width 


Distal  calibre  (IQrnrn  from  extremity)  = 
Distal  width  x distal  thickness 


Robu  st  icity  index  = 
Mesial  width  ^ 

Maximal  length 


Afirjlf* 


Caviiydiamei  er 


Thickness  of  The  base 


Q £ Width  d the  t?gw 


Internal  diameter 

< > 


Entennal  diameter 


irUerorbilal 

breadth. 


^reatkh  orhrainca^c 


iJLorrHta]  brcacth 


Nasa;  width 


Greatest  length  of  skull 


/.  c/e  ft'Cxn  z — ^ x - 3 3 ’§ 

*3 


Curvas  para  medicion  de  raspadores.  y ejemplo  de  su  aplicacion  a un  raspador  de  frente 
mixto  o complejo. 


Fig.  378 


What  is  Shape? 


• a shape  can  be  understood 
as  any  “single”, 

“distinct”,  “whole”  or 
“united”  visual  entity. 


What  is  Shape? 


Properties  of  a configuration  of  points  which  are 
not  altered  for  effects  of  size,  position  and 
orientation,  or  by  translation,  rotation  and  scaling 


^ WILEY 

Statistical 
_ Shape 
Analysis 


I 


Las  L Drvden  ,ni>  Kwn  V.  Mamma 


Describiendo  el  Perfil 


Fig.  6.2 


£chantillonnage  de  profit 


SJ.  Lycett  ef  al.  / Journal  of  Archaeological  Science  _o  (2005)  J— J5 


Fig.  6.  Bilateral  measurements  taken  from  the  length  axis  of  the  base.  A total 
of  26  measurements  were  taken  at  percentage  points  left  (A)  and  right  (B)  of 
this  line. 


Fig.  7.  Using  setsquares  to  accurately  align  and  position  uprights  along  iht 
scales  of  the  base. 


Fig.  I.  The  steps  leading  from  a scanned  drawing  to  the  curvature  function,  a.  A scanned  drawing  of  a bowl.  b.  The  pixelized  profile,  c.  Enlarged 
detail  of  the  pixelized  profile  and  the  interpolated  curve,  d.  The  curvature  function. 


Fig.  1*  Single  profiles  of  two  different  bowls  from  Tell  Dor.  Israel  (a.c 
and  the  overlap  of  six  profiles  of  the  same  bowls  (b.d).  respectively 


Relational  Indexes 


Circularity:  the  degree  of  circularity  of  a texel.  That  is  how  much  this  texel  is 
similar  to  a circle.  Where  1 is  a perfect  circle  and  0.492  is  an  isosceles  triangle. 
This  shape  is  expressed  by: 

4 7ZS  s:  texel  area 

— p:  texel  perimeter 

P~ 

Quadrature : the  degree  of  quadrature  of  a texel.  where  1 is  a square  and  0.B00 
an  isosceles  triangle.  This  shape  is  expressed  by: 

P 


2 M 


Irregularity:  measurement  of  the  irregularity  of  a texel.  done  with  its  perimeter 
and  the  perimeter  of  the  surrounding  circle.  The  minimum  irregularity  is  a circle, 
corresponding  at  the  value  1.  A square  is  the  maximum  irregularity  with  a value 
of  1.402.  This  shape  is  expressed  by: 

Pc_ 

P 

Elongatioru  as  the  degree  of  ellipticity  of  a texel.  Where  a circle  and  a square  are 
the  less  elliptic  shape.  [This  shape  is  expressed  by: 


D 

d 


D:  maximum  diameter  within  a.  texel 
d:  minimum  diameter  p Crp Clldl C ular  at  D 


bw 


Image J . http : //rsb . inf o . nih . go v/ij  / 


L Image  J 


File  Edit  Image  Process  Analyze  Plugins  Window  Help 

[3 olgjcjSNKI-fKlAl  KlnW 


Location  = (1 44,85),  value=D,66,200 


1 67x1 58  pixels;  8-bit  grayscale;  :j 


212x151  pixels;  RGB;  124K 


- >3  - 


- 


-Inj  *1 

121x118  pixels;  8-bit  gr; 


-Ini  *1 


3E  Jsi-sJ 

102x194  pixels;  8-bil 
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Results 


File  Edit 


Area 

Mean 

Major 

Minor 

Angle 

1 

425 

195.95 

23.02 

19.31 

71.22 

2 

426 

201.84 

31.33 

17.31 

17.59 

3 

676 

198.99 

35.72 

24.10 

166.25 

4 

361 

197.21 

23.70 

19.39 

172.83 

5 

610 

189.72 

46.20 

16.81 

64.39 

6 

641 

192.62 

39.75 

20.53 

122.64 

Image J . http : //r sb . inf o . nih . go v/ij  / 

Fiji  is  Just  ImageJ. 
http://fiji.se/Fiji 


0 0 Set  Measurements 

[Vj  Area 

'V*  Mean  Cray  Value 

_ Standard  Deviation 

'12  Modal  Cray  Value 

_ Min  & Max  Cray  Value 

_ Centroid 

'_J  Center  of  Mass 

Perimeter 

_ Bounding  Rectangle 

_ Fit  Ellipse 

_ Circularity 

□ Feret's  Diameter 

_ Integrated  Density 

_ Median 

Skewness 
5?  Area  Fraction 

_ Kurtosis 

_ Limit  to  Threshold 
_ Invert  Y Coordinates 

Display  Label 

Redirect  To:  None  i 1 

Decimal  Places  (0-9):  1 

( Cancel  OK 

a\ 

Consider  particles  having: 

Particle  size: 

0,99999999 

Aspect  ratio  (>1): 

1.0,50 

Form  Factor  (0-1): 

0.0, 1.0 

Compactness  (0-1): 

0.0, 1.0 

Roundness  (0-1): 

0.0.1.0 

Angle  (0-130): 

0.0,180 

'V*  Clear  particles 

Wj  Include  holes 

Calculate 

Single  particle  details  * 1 

Decimal  places: 

2 

'Vi  Display  orientation  plot 

Weigh  orientation 

Interval  width  (0-40  degrees): 

10 

2J  Include  edge  particles 

( Cancel  OK  ) 
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Measuring 

Shape 


An  ordered  n-tuple  of  measurements  completely  characterizes  a shape 
without  redundancy  if  (a)  there  is  a set  of  drawing  rules  that  permits 
reconstruction  of  the  shape  outline  using  only  this  ordered  n-tuple  of 
measures,  and  (b)  there  is  no  ordered  /(-tuple  of  measures,  k<n,  such  that  the 
shape  outline  can  be  reconstructed  from  the  ordered  /(-tuple  (Read  2007). 


Morphometries 


ikklirthL  ftlLllH  111! 


Morphometries  for 

Non  morph  ometricia  ns 


t.  Spring 


Landmark-based  geometric 
morphometries : summarises  shape  in 
terms  of  landmark  configuration 
(discrete  anatomical  loci,  described  by 
two-/three-dimensional  Cartesian  co- 
ordinates); 

Outline-based  geometric 
morphometries:  summarises  the  shape 
of  open  or  closed  curves,  typically 
without  fixed  landmarks.  Analyses 
include  Fourier  (Elliptical  Fourier,  Fast 
Fourier  etc.)  and  Eigenshape  (and 
Extended  Eigenshape)  analyses 


ArchaeoMorph 

www.archaeologicalmorphometrics.com 

L J 
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Morphometries,  the  practical  side 


tpsUtil  and  tpsDig  programs  from  Rohlf ’s  morphometries  site 

http://life.bio.sunysb.edu/morph/ 

Point  Picker  plugin  for  ImageJ 

http://bigwww.epfl.ch/thevenaz/pointpicker/ 

MorphJ.  Manchester  University 

http://www.flywings.org.uk/morphoj_page.htm 
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LOHSE,  E.S.,  SCHOU,C.,  SCHLADER,R.,  SAMMONS, D,  2004,  “Automated  Classification  of  Stone 
projectile  Points  in  a Neural  Network”.  In  Enter  the  Past.  The  e-way  into  the  four  dimensions  of  culture 
heritage.  Edited  by  Magistrat  der  Stadt  Wien-Referat  Kulturelles  Erbe-Stadtarhchaologie  Wien.  Oxford, 
ArcheoPress  (B  AR  Int.  Series,  S1227),  pp.  431-437). 


Tip  Landmark 


Fig.  3.  Image  of  Clovis  point  showing  approximate  location 
where  1 1 interlandmark  characters  are  measured  and  the  location 
of  the  three  landmarks.  Character  initials:  EL,  edge  boundary 
length;  TB,  tip  landmark  to  base  landmark:  TW,  width  of  tip  to 
base  length  to  maximum  inflection  position;  BL,  blade  length; 
MW,  maximum  width;  BB.  base  boundary  length;  LB.  linear 
measure  of  base;  ML,  mid  line  length;  OL,  overall  length;  BW, 
basal  width  across  first  third  of  point;  LT,  length  from  base  to  1/3 
along  opposite  edge. 


End  Point 


Fig.  8.  S- shaped  vessel:  profile  segmental  ion  scheme. 


X 


*- 


• MI,  local  minimum:  point  of  vertical  tangcncy;  point 
where  the  jc- value  is  smaller  than  in  the  surrounding 
area  of  the  curve. 

• MA,  local  maximum:  point  of  vertical  tangcncy;  point 
where  the  .v-value  is  bigger  than  in  the  surrounding  area 
of  the  curve;  they -value  refers  to  the  height  of  the  object 
(e.g.  MAO?)). 

• CP,  corner  point:  point  where  the  curve  changes  its 
direction  substantially. 

• BP,  base  point:  outermost  point,  w'hcre  the  profile  line 
touches  the  base  plane. 

• RP,  point  of  the  axis  of  rotation:  point  wrhcrc  the  profile 
line  touches  the  axis  of  rotation. 

• OP,  end  point:  point  where  the  profile  line  touches  the 
axis  of  rotation;  applied  to  incomplete  profiles. 


3 dimensionality  of  Physical 

Space 


z 
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SHAPE  DECOMPOSITION 

Geon  Theory  (Biederman  1987) 


ComuKtcd- 

Ori^iniiL 


Disconnected- 

Ori^maJ 


Curnitcted-  Disconnected* 

Scrambled  Scrambled 


V 

StraigM  Edge 
EirdifEi  Axis 
Cora&nr 


Aith 


3tai£iv!  Edge 
Qirted  Axis 
Cora  twit 


Wedge 


Straisii;  Edge 

Su  aigti’  A-j's 
Enpiwipd 


Ooce 

A 

Curved 
SUdiHh:  .Axi  s 
Ewstdpd 


FVramifl 


Sire  ight  Tiitc- 
Sunlit  AJiifl 

EipiivM 


Cyiiii'lei 


Cnrofl  Edge 
Sb^lil  taia 
Expand 


Curved  Ed&e 
Sb  - j i s"  El  Axis 
Conutasi 


Middle 


Curved  E3«$e 
Curved  Ass 
Con?t3firt 


Barrel 

O 

Curved  Edge 
;3u  ii>'ln.  .-jas 

S^j  &.Ccn‘ 


Expanded 

ILdnJh 


Curved  Edge 
CuiVgd  A Kid 


3D  Morphometries 


Center  for  Regional  Heritage  Research 
Robert  Z.  Selden  J r¥  Ph.D*  RPA 

Calderon  Collection  - 11-005  jO 
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*U9fil 
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Photogrametry 


oriented  images 


Photogrametry 


Perspective 


faces:  1000000  vertices:  500061 


iijx  ^ ^ a g a §& 
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Photogrametry 
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3D  Scanning 


3D  Scanning  / 


Target  reflects 

linht  tn  hhoto  Cell 


Photocell 


Oscillating 
mirror 


Laser 


CJtteehhg  Lara 

X - ^ijJ  Scanner 


FiMd  Mmyet 


V ^ *,xi*  S^rw^i' 


J-ihkJ  M rrof 


awirror, 

V 
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3D  Scanning 
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3D  Scanning 


BREUKMAN  OTO 
Institucio  Mila  i Fontanals 
(CSIC) 


Minolta  Vivid 
UAB-CVC 
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3D  Scanning 


0 NextEngjne 


HEME  r.KJDUC 


AsCNfiflBrtne 

The  #1  Selling  30  Scanner. 

TI4  30  &«ifwr  objMt  in  hill  »tor  with 

m utti- LnM-r  S3*raojson  B'oeirttTo^flhioGhrolc-gyhas  ?m.3da  rubs 
World  Si  rriEri!  poptjlnj  2f>  Scjnnnr.  win  iboLKnndc  Wf^rs.  fi  av-nr  £0 
COunUtM-. 

tn  fine  bow  ovosythng  you  rmed  Da  cUgtiza  3D  models,  InaOudnig 
ScanSludio  HQ -bjlt^re.  Ewporlilo  STL  OfiJ  YftML  XT^anc 
otter  ronrisas. 

Otfput  3D  scan  iicxlehs  l&pco  j ji  assign  wlforare  l-Ui  golltfWofk-t. 
3d sfifla*  ZBrusi  Rhino 

SeanStudm  CAD  Toot*  Id  ^§eklymjtosrffae&  ftes  r>j  fophfWoikf- 
DC  hut:l  r Ftps  Prr.1  iri  Oimon^ign  W Syplflm*  iCgip 

Objol  ;v:d  olhrrr  HD  jsrmlr"- 


n yh  ptrfciiii  allot  dl  UlCHh  Nit  pli'^T 


HU 
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3D  Scanning  / Photogrametry 


DavidScan 


2995$  AMAZON.com 
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3D  Scanning 


SCANIFY 


€1199.00 


price  excluding  tax 


Buy  now.. 


3D  scsnnirtg 
made  easy. 


You’ve 


FLrti.%0 
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Louer 


3D  Scanning 

1FARO  FOCUS  3DJ 


Dia 

Semana 

15  dfas 

Mes 

Laser,  trfpode  y 
software  Scene 

250 

860 

1500 

2550 

8 Esferas 

10 

30 

50 

80 

- Clase  practica  inicial:  110  € (2  horas) 

- Seguro:  incremento  4% 

- Tarifas  validas  durante  2015 

- IVA  no  incluido 


Caracteristicas: 

• 5kg  de  peso 

• Velocidad  de  976.000  puntos/s 

• Precision  ± 2mm 

• Camara  integrada  de  70  megapixel 


XYZ 
113 
436 
7 6 14 


XI  X2  X3  X4 
Y1 
Y2 
Y3 


127/  81 


Formas  tridimensionales  = Interpolation 


Fignra  2.20  Malha  paligoin^L  com  10.000  (es^.)  e 1.250  (drtj  poligonos  f^loitiiLb&OVb,  c] 


Fignra.  2_iQ  FJor^arir^g^plicada  3 malha  poligoiLs]  coni  10.000  (etq  ) e 1.250  (drt.)  poligono*  [5,lcdtinno0  7b.  c]. 


Fignra  2.J1  Smooth  shading  aplic;i{b  a rr  a1:«i^  poligazui  com  10.000  (etq.)  e 1.250  (drt.)  polisonoi  p^aianhoOTb.  c]. 
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3D  File  Format 


facet  normal  ni  rij  nk 
outer  loop 

vertex  vl  vlv  vl 

x y z 

vertex  v2v  v2„  v27  vertex  v3v  v3v 

x y z x y 

v3z  endloop 
endfacet 


3D  Shape  Descriptors 


3DVIA  MP  - Test  / Project  Editor 


d & y 


^7  Presets 

lA  2D 

sV3D 

Cameras 

CDSDynamics 

B ^ Characters 

El/y  Input  Devices 

9 Lights 

* i Misc 

S*"  Particle  Systems 

El5#  Physics 

Shaders 

^7  Project 

li 

E1&  UpdateTest 

^7  User 

la 

< a ^ na 


'W  Physical  Collision  Group 

Add  an  object  group  to  define  new  collision  rules  (a  physical  worl 

■ra 

Physical  World 

Add  a physical  world  which  allows  interactions  between  physical  ■ 

W Physical  Zone 

Add  a fixed  zone,  other  objects  can  generate  collisions  events  wi 

Physical  Rigid  Object 

Add  a physical  behavior  tc 
Sy  world  is  required) 


s 3D  object  (must  be  dropped  on  this  3D  object,  a physical 


Project  Editor  fr-ymmi  , 

A u 


I Cfg  | Comps [ Present  | Status  | 


E)  Permanent  Stages 
□ •a  DefaultStage 

Dependencies 
Resources 
El  Templates 
B Actors 

Main  Camera 
Hi  Viewport 
Physical  World 
E ^3  Model 

meshl 

v Omni  Light 
Dynamic  Stages 


Si  Ready 


/ 

/ 


T 


2) 

V 


ESBlilii  C Selection  : meshl(vkPhysicsCDSDynRigidObject) 
vkPhysicsCDSDynRigidObject 


Component  D 

Active  At  Present 

Active 

E Representation 

Graphical  Object 
Move  Graphical  Object 
Representation  Scale 
Execution  Visible 
B Object  Specifications 
El  Shape 

Bounding  Shape 
El  Bounding  Position 
El  Bounding  Orientation 
Radius 

Shape  Change  Detection 

Density 

Mass 

El  Gravity  Center 
El  Inertia  Momentum 
El  Inertia  Orientation 
Friction 


El 
SI 

El 

~ 1.000 

□ 

Type  : * vkPhysics CDS DynSp here 
Type  : - *null 
Edit...  0.000.0.000.0.000 
Edit...  0.000.0.000.0.000 


- 0.044 


..  0.000 
3l.000 

Edit...  0.000.0.000.0.000 
Edit...  1.000,1.000.1.000 
Edit...  0.000.0.000.0.000.1.000 
Q 0.400 
~ n snn 
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3D  Shape  Descriptors 


Sphericity  f =(0/0^ 


Roundness  R = (2r:/Dj) 


Shape  factor 


Circularity  Factor 

Circle  of  equal 
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CUBENESS 


^Cubcncsft  measure  C^=31J.d(5)  cannot  actually  reach  zero 


M 


but  the  computed  value  i&  0.0000  due  to  numerical  error. 


4Cubetiess  measure  Cpm  ] s.o( S ) cannot  actually  reach  zero  but  the  computed 
value  is  0,000  due  to  numerical  error. 
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Texture 
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Texture 


Definition  of  those  attributes  of  a surface 
that  have  visual  or  tactile  variation  and 
contribute  to  distinguise  those  surface  from 

others  interacting  with  it 


Any  surface  shows  variation  in  its  properties  like:  albedo, 
color,  uniformity,  density,  roughness,  regularity, 
linearity,  direccionality,  brightness,  defomation, 
reflectivity,  opacity,  transparency, 


Texture 


COLOUR: 


The  quality  of  an  object  or  substance  with 
respect  to  light  reflected  by  the  object, 
usually  determined  visually  by  measurement 
of  hue,  saturation,  and  brightness  of  the 
reflected  light;  saturation  or  chroma;  hue. 


Texture 


MAGENTA 


AZUL 


Circulo  cromatico  RGB 
(LUMfNICO) 


AMARILLO 


Colores  primarios: 
Rojo,  Verde,  Azul. 
Colores  secundarios: 
Cian,  Magenta,  Amarillo. 


VERDE 


o •: 

c£-  JD 

••• 


RANJADO 


OLACEi 


AMARILLO 


VERDE 


MAGENTA 


Circulo  cromatico  CMY 
(PIGMENTARIO) 
Colores  primarios: 

Cian,  Magenta,  Amarillo. 
Colores  secundarios: 
Violaceo,  Anaranjado,  Verde. 


Texture 


Hue 


Saturation 


Brightness 


Texture 


Value 


Texture 


COLOR:  Colorimeter 


Como  funciona 
un  colorfmetro 


Fillies de  folodeleclore-s  Los  valorem 

colores  que  oonvierten  se  process  n 

primaries  la  I uz a va  lores  a va lores  Lab 

aditivos  numericos  conforms  a un 

lluminante  estandar 


Muestra 

de  luz  (par  emisidr 
o refleitidnl' 


Resultado: 
Valorem  Lab- 

del  color  recibido 


TACTILE  APPEARANCE 


High  Roughness 


Moderate 

Floughr&s& 


>750 


550-750 


Low 

Roughness 

350-550 


Smooth 

200—350 


Polished 

<200 


Flaw 

( unspecified) 


r 


Nomina  I surface 
Nonna!  section 


f y // _////;  ///  /////  // / // ////  / 

J4C^rLrlfcC/  / //  ///  //  // / //  //y  * 

i' 

Total  profits 

^mnr  in 
gctsmelric  forml 

Wa  vines®  pi  n ti  le 
(roughness  bcighis 
a nconatcdli 

Roughness  profile 
iwavincss  heights 
atienuased) 

Coarseness, 

roughness, 

smoothness, 

polish, 

burnish, 

bumpiness, 

waviness 
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Wear  traces  as  texture 


In  use  wear  analysis,  we  are  interested  in  the  nature  of  the  alterations 
suffered  on  the  lithic  tool  surface.  In  this  way,  our  aim  is  to  make 
inferences  about  the  working  processes  developed  by  people  that  used 

these  tools 


Original 

Surface 

Altered 

Surface 


r:;HV5S- 


Decoration  as  Texture 
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Surface-microtopography: 

Coarseness,  roughness,  smoothness,  polish,  burnish,  bumpiness,  waviness 


VISUAL  TEXTURE 
is  a 

GEOMETRICAL  MODEL  OF  LUMINANCE 

we  recognise  it  as 
IMAGES 


33333*45 5 AG666&G 
4333a45G66£&fi6B6 
4 3 3 3 4 5 7 6 3B  7 S A 
3322  & 7 N (}£  0G5  55  &i 
33Z257BB9B7SA555 

2 1 224$  £9967  &4  4 &5 

2 I 1 235649# 7 54444 
1101476903754444 
l 1 0 n 63  fl*H  ^ g 5 & TJ  l 
1Q0T3579QB746S55 
■;1  Q Q 1 2 57  ft#  7 7fi  5 56  & 
1 T 1 237777655555 
1 11  13S6?775555ia 
1 1 1 22  J 5 6tl  55  ± 55  55 
2_2  2 3 2 4 4 SB  55  5 fl  % J S 


1)  IMAGE  DECOMPOSITION 

• We  look  for  luminance  intervals, 
to  select  those  pixels  in  the  image 
(through  a density  slice),  which 
coincide  with  observed  texture 
discontinuities 

• These  three  intervals  are  those 
which  allow  a better  description 
of  texture  patterns 


lithics 


0-80 


0-120 


160-200 


Original  Image 


Segmentation  “Texture  Components’ 

1 


DESCRIBING  TEXTURE 

The  features  that  we  take  in  account  are: 


By  Shape: 

- Elongation 

- Circularity 

- Quadrature 

- Thinness 

- Ratio  Compactnes/Thinness 

- Compactness,  mesured  through  two 
equations 

- Irregularity 

- Rectangularity,  mesured  through  two 
equations 

- Ratio  Perimeter/Elongation 

- Feret  diameter 

- Minimum  rectangularity 


By  Composition: 

- Mean,  mean  of  luminance 

- SD,  standard  deviation  of  luminance 

- Mode,  mode  of  luminance 

- Min,  minimum  luminance  value 

By  Size: 

-Area 

-Major  axis 

-Major  axis  perpendicular  to  the 

major  axis 

-Perimeter 


Surface-microtopography: 

Coarseness,  roughness,  smoothness,  polish,  burnish,  bumpiness,  waviness 


Elevation 

(m) 


4 0 
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CON-FOCAL  MICROSCOPY  AT  DIFFERENT  SCALES 


Light  source 


multiple  X,Y 
"scans  per  Z plane 


/<*  o o o/ 

AJOOO/ 

/o  o o 
/o  ooo 

y 

1 ZOOOI 

/booo 

7 

3 Z 

movements 


Confocat  Microscopy 


Light  source 


Z 


z 


All  X,Y  data  captured 
in  single  acquisition 
per  Z plane 


3 Z 

movements 


WLI  3D  Microscopy 


of  Celt 
Layer 
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Leica  DCM8  METROLOGY  MICROSCOPE  OF  3d  SURFACES  de  superficies  3D 
combining  confocal  microscopy  and  interferometry 
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Surface-microtopography: 


Coarseness,  roughness,  smoothness,  polish,  burnish,  bumpiness,  waviness 


A modern  range  scanner 
captures  surface  data 
points  less  than  50 
microns  (0,05  mm),  apart 
from  producing  high- 
density  triangular  meshes 
with  an  average  resolution 
of  over  1000  points  per 
cm2. 


iz j. ijm) 
h?.  i* 
IH.&t 

int .34 

33.33 
9% 00 
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Surface-microtopography: 

Coarseness,  roughness,  smoothness,  polish,  burnish,  bumpiness,  waviness 


Mechanical  profilometer: 
Nanometric  resolution 
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Surface-microtopography: 

Coarseness,  roughness,  smoothness,  polish,  burnish,  bumpiness,  waviness 


Profilometer 

Resolution:  pm  / 0,006  pmlOO  jim  / 0,002 
pm25  pm  / 0,0004  pm 


127/109 


Surface-microtopography: 

Coarseness,  roughness,  smoothness,  polish,  burnish,  bumpiness,  waviness 


Alta  precision  en  el  tipo  1.7pm 
Crysta  Apex  S es  una  MMC  CNC  que 
garantiza  un  error  maximo  admisible  de 
medicion  MPEE  = (1.7  + 3 L/1000)  pm  [Serie 
500/700/900]. 

Sistema  de  Compensacion  de  Temperatura 
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Surface-microtopography: 

Coarseness,  roughness,  smoothness,  polish,  burnish,  bumpiness,  waviness 


Waviness  measurement 


Gaussian  Curvature 
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3D  view  of  surface  texture 

with  form  removed  pseudo-colour  view 


Ah  3D  Digital  surface  texture  patterns  description 


original  subtracted  surface  texture  after 

surface  form  form  subtraction 


sample  Ast_s09 


sample  Ast_s01 


sample  Ast_s09 


sample  Ast_s01 


Ah  ISO  25178  series 


ISO  (2012).  ISO  25178-2:2012  - Geometrical  Product  Specifications  (GPS)  - Surface  texture:  Areal  - 
Part  2:  Terms , definitions  and  surface  texture  parameters.  International  Standard  Organization  (ISO). 


The  ISO  25178  series  define  more  than  40  areal  parameters,  grouped  in: 

• Height:  statistical  distribution  of  height  values  along  the  z axis. 

• Spatial:  spatial  periodicity  of  the  data,  specifically  its  direction. 

• Hybrid:  spatial  form  of  the  data,  i.e.,  amplitude  and  spatial  information. 

• Functional  volume:  surface  bearing  area  ratio  curve  (Abbott-  Firestone  curve). 

• Features:  selected  features  are  identified  by  segmentation. 


Calculations  are  made  upon  the  entire  surface 


AND  NOT  upon  averaging  estimation  calculations  derived  from  2D  profilometric  methods  and 


ijh’ 

Computational 

Surface 

and  Roundness 
Metrology 


Richard  Leach  Editor 


Characterisation 
of  Areal  Surface 
Texture 


(ASME  : 


L-iL-jr-d  E-dUnn 


Surface 
Topography 


k j siout  & i Biuni 


Pinion  | 
Press 


Springer 


Sa  and  Sq 


Sa  and  Sq  are  the  Average  Roughness  and  Root  Mean  Square  Roughness  are  evaluated  over  the 
complete  3D  surface  respectively.  Mathematically,  Sa  and  Sq  are  evaluated  as  follows: 


\\\Z(x,y)\dxdy 


Plateau-like  surface  Sa  = 16.03  nru  Sq=  25.4  nm 


Sq  ~ ^Jjj(Z(y\y)f  drtlr 


Surface  with  Peaks  Sa  = 16.03  nm  Sq=  25.4  nm 


127/ 


114 


Ssk  (Skewness)  and  Sku  (Kurtosis) 


Ssk  and  Sku  are  the  Skewness  and  Kurtosis  of  the  3D  surface  texture  respectively  Figuratively,  a histogram  of  the  heights  of 
oil  measured  points  is  established  end  the  symmetry  and  deviation  from  an  ideal  Normal  (i.e.  bell  curve)  distribution  is 
represented  by  Ssk  and  Sku  Mathematically,  the  Ssk  and  Sku  are  evaluated  as  follows: 


s 


jj  (Z(;vs  y)fcbcdy 


Surface  with  multiple  peaks  Ssk  = 3.20  Sku  = 18.71 


Periodic  Texture  Ssk  = 0.16  Sku=  1.63 
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Sp  (Max  Peak  Height),  Sv  Max  Valley  Depth)  and  Sz  (Max  Height  of  Surface) 


Sp,  Sv;  and  Sz  are  parameters  evaluated  from  the  absolute  highest  and  lowest  points  found  on  the  surface  Sp,  the  Maximum 
Peak  Height,  is  the  height  of  the  highest  point , Sv,  the  Maximum  Valley  Depth,  is  the  depth  of  the  lowest  point  (expressed  as  a 
negative  number)  and  Sz  the  Maximum  Height  of  the  Surface),  is  found  from  Sz  = Sp  - Sv. 

Note:  earlier  standards  referred  to  Rz  as  a average  of  the  10  highest  to  10  Lowest  Points  and  other  variations.  The  ISO 
community  agreed  for  the  newer  standard,  ISO  25178-2  to  establish  Sz  as  strictly  the  peak  to  valley  height  over  a a real 
measurement. 


A surface  used  in  the  printing  industry  characterized  by  deep  valley 
structures  with  Sv  being  - -15pm 


A.  polymer  surface  prepared  with  asperities  as  measured  by 
Sp  being  ~ 0.90  pm 


ACF  (Autocorrelation  Function) 


Rpm,  Rvmi  and  Rz  (Average  Max  Peak  Height,  Valley  Depth  and  Height  of  Surface) 

The  Rpm,  Rvm,  and  Rz  parameters  are  evaluated  from  an  average  of  the  heights  and  depths  of  a number  of  extreme  peaks 
and  valleys.  Rpm  the  Average  Maximum  Peak  Height,  is  found  by  averaging  the  heights  of  the  ten  (10)  highest  peaks  found 
over  the  complete  3D  image  Rvm,  the  Average  Maximum  Valley  Depth,  is  found  by  averaging  the  depths  of  the  ten  (10) 
lowest  valleys  found  over  the  complete  3D  image  Rz,  the  Average  Maximum  Height  of  the  Surface,  is  found  from  Rpm-Rvm 
Note  that  in  determining  the  peaks  and  valleys,  the  analysis  software  eliminates  a grid  of  1 1 x 1 1 pixels  around  a given 
peak/valley  before  searching  for  the  next  peak/vailey,  thus  assuring  that  significantly  separated  peaks/valleys  are  found. 


imm 


^bih  ilium 


* ^.uu 


127/ 


117 


a 


b 


c 


d 


Actin  orientation  histogram 


Actin 
single  cell 


Actin 

o ri  ertation  colou  r code  d C rop 


Primary  oriental  ion 
Secondary  orientation 
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SURFACE  METROLOGY 


http : //www.michmet . com/index,  html 

Free  Software:  http://digitalmetrology.com/free-software/ 

The  freeware  version  of  SigmaSurf  allows  users  to  import 
data  (through  an  ASCII  file  format)  and  study  the  effects  of 
various  filters  while  viewing  the  Primary,  Waviness  and 
Roughness  profiles  at  various  scales. 
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Using  a 3D  scan  for  measuring  surface- 
microstructure.  Neolithic  wood  tools 


j 10  mm 


Fragmentation  and  restoration  techniques  (surface 
finishing):  holes  and  noise  data. 


Some  difficulties  during  3D  data 
capture: 


Surface  area  with  distinct  characteristics  - restoriation 
product  (A),  wood  hardened  with  fire  (B)  and  natural  wood 
(C):  holes  and  noise  data. 


Vera  Moitinho  de  Almeida,  Juan  Anton  Barcelo  | 1st  International  Conference  on  Best  Practices  in  Cultural  Heritage:  Archaeology  | Menorca,  09-13. April. 2012 


12  /25 
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Laboratori  d’Arqueologia  Quantitativa  ^HABTTQlfl&tjQn 


urB 

Universitat  Autbnoma 
de  Barcelona 
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Average 

Roughness 

Sa  = (1 /At)  j'  1 lZ(x,y)\  dxdy 

Root  Mean 

Square 

Roughness 

Sq  = [(1  /Ae  jj  jo  Z*Or,yJ  dxdyj  * 

Skewness 

M’^zCCz>’^ 

Kurtosis 

zV,-!"“y 

Texture 

Aspect 

Ratio 

, length  of  - fastest  - decay  - AACV  -in  - any  direction 

Length  of-  slowest  - decay  AAdV  -in-  any  direction 

Cubeness  Ca{S)  - "(S)  ~i4(5)/e 

* ' n-(yn(S))» 


Direction  Major -direction -of  - Lay -Derived- From -APSD 

Surface 


Ah  3D  Digital  surface  texture  sampling 


Experimental  surfaces 

10  samples  (50*50  mm)  [E] 


Experimental  stela  (half  sculpted) 

6 samples  (50*50  mm)  [Est] 


Archaeological  stela  (SMB/08  e-17/5/958,  fragment  a) 

27  samples  (50*50  mm)  [Ast] 


3D  structured  light  scanner 
(SmartSCAN3D  Duo  System, 
Breukmann) 

FOV:  450  mm  stereo 
Resolution:  280  pm 

(according  to  manufacturer) 


And  now,  what  we  can  do  with  all  that? 


WOULD  YOU  RELAX 
ALREADY  T?  YCH/TRfc 

. OM  vacation:  „ 


The  usual  archaeological  answer 
typology  and  classification 


But  Typology  is  not  Classification 

Archaeological  “types”  do  not  give  “answers  nor  explanations 


Was  it  used  as  a bottle? 


Was  it  used  as  an  arrow  point? 

YES 


Learning  Stage 


Partitioning  Stage 


THE  TASK: 

to  find  the  common  structure 
in  a given  perceptual 
sequence 


ASSUMPTION: 

the  structure  that  is  common  across 
many  individual  instances  of  the  same 
cause-effect  relationship  must  be 
definitive  of  that  group 


From  resemblance  to  similarity 

Studying  “resemblance” 

Two  entities  are  similar  because  they  have  many 
properties  in  common.  According  to  this  view: 

• similarity  between  two  entities  increases  as 
a function  of  the  number  of  properties  they 
share 

• properties  can  be  treated  as  independent 
and  additive 

• the  properties  determining  similarity  are  all 
roughly  the  same  level  of  abstractness; 

• these  similarities  are  sufficient  to  describe 
a conceptual  structure:  a concept  would  be 
then  equivalent  to  a list  of  the  properties 
shared  by  most  of  its  instances. 


Understanding  “Similarity” 

Studying  “resemblance” 

The  very  idea  of  similarity  is  insidious.  First,  we 
must  recognize  that  similarity  is  relative  and 
variable.  That  means  that  the  degree  of  similarity 
between  two  entities  must  always  be  determined 
relative  to  a particular  domain.  Things  are  similar 
in  color  or  shape,  or  in  any  other  domain.  There  is 
nothing  like  overall  similarity  that  can  be 
universally  measured,  but  we  always  have  to  say 
in  what  respects  two  things  are  similar.  Similarity 
judgments  will  thus  crucially  depend  on  the 
context  in  which  they  occur. 


Studying  “resemblance” 


Measuring  “Similarity” 


• Dice  (Sorensen)  coefficient  for  absence-presence  (coded  as  0 or  positive 
numbers).  Puts  more  weight  on  joint  occurences  than  on  mismatches. 

When  comparing  two  columns  (associations),  a match  is  counted  for  all  taxa 
with  presences  in  both  columns.  Using  ’M’  for  the  number  of  matches  and 
’N*  for  the  the  total  number  of  taxa  with  presences  in  just  one  column,  we 
have 

Dice  similarity  = 2 M /(2 M + N) 

• Jaccard  similarity  for  absence-presence  data:  M /(M  + N) 

• The  Simpson  index  is  defined  as  M /Nmur.  where  Nrriin  is  the  smaller  of 
the  numbers  of  presences  in  the  two  associations.  This  index  treats  two  as- 
sociations as  identical  if  one  is  a subset  of  the  other,  making  it  useful  for 
fragmentary  data. 

• Kulczynski  similarity  for  presence-absence  data: 

M/(M  + Ni)+  M/(M  + JV2) 

2 

• Ochiai  similarity  for  presence-absence  data  (binary  form  of  the  cosine): 

+ Ni)][M/(M  + iV2)] 

• Bray-Curtis  measure  for  abundance  data. 


Bray  — Curtis ^ = 


E?=i  kij  - xik  I 
Ei=l  {Xij+Xik) 


• Cosine  distance  for  abundance  data  - one  minus  the  irmer  product  of  abun- 
dances each  normalised  to  unit  norm. 


Studying  “resemblance” 

• Measuring  distance  as  an  alternative  to 
“similarity” 


From  Clustering  To  Classification 


• The  Clustering  Principle 


The  Clustering  Principle: 

Internal 


Similarity 


higher  than 

External 


Similarity 


Law  of  Similarity: 

Items  that  are  similar  tend  to  be  grouped  together 

In  the  image  above,  most  people  see  vertical 
columns  of  circles  and  squares. 


Clustering 


It  is  the  process  of  grouping 
input  samples  in  similarity 
classes. 

Thus,  one  would  like  to  group 
observations  so  as  to 
minimize  intra-group 
distances  while  maximizing 
inter-cluster  distances,  subject 
to  the  constraints  on  the 
number  of  clusters  that  can  be 
formed. 


5 ® 

LuoQhezza 

0 15  20  25  30  35  40  45  50  55  60 


Linear  Discrimination 


;■  f ( D.  '-^rPr'e  \ 95$  '2r. 


More  generally 

f(x)  = H WiXi  - 0 

= |w||x|  COS  0-0 


for  cos  0 

-90  to  +90— ►class  1 
else— ►class  2 
ie.  for  +1  or  -1 
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Non-Linear  Discrimination 





From  Clustering  to  Classification 


• What  is  a “Classification”? 


Classification  is  a form  of  categorization 
where  the  task  is  to  take  the  descriptive 
attributes  of  an  observation  (or  set  of 
observations)  and  from  this  to  identify  the 
observation  within  a different 
phenomenological  domain.  Hence,  the  task  of 
the  classifier  is  somehow  to  partition  feature 
space  into  disjoint  regions  that  each  represents 
a particular  class,  cluster,  or  pattern. 


From  Clustering  to  Classification 


• What  is  a 


“Classification”? 


The  goal  in  a classification  problem  is  to 
develop  an  algorithm  which  will  assign  any 
artifact,  represented  by  a vector  x,  to  one  of  c 
classes  (chronology,  function,  origin,  etc).  The 
problem  is  to  find  the  best  mapping  from  the 
input  patterns  (descriptive  features)  to  the 
desired  response  (classes). 


From  Clustering  to  Classification 


• What  is  a “Classification”? 


The  purpose  of  the  classification  problem  is  to 
estimate  the  probability  of  membership  of  the 
case  in  each  class.  The  objective  is  to  build  a 
model  with  significant  predictive  power.  It  is 
not  enough  just  to  find  which  relationships  are 
statistically  significant.  That  explains  why 
classification  and  prediction  are  frequently 
interrelated.  A prediction  of  an  historical  event 
is  equivalent  to  a classification  within  a given 
set  of  events. 


COMPUTER  VISION 
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3-  Scene  Analysis  Unit 


Feedback 


Stages  of  Computer  ’vision,  ©2005  F.  Keynoua 


4-  Visual  ftrception 


Input  = (10  42  6 .1) 

ICONOGRAPHIC 
COMPLEXITY  = 10 

RELEVANCE  OF  HUMAN 
FIGURE  = 4 

RELEVANCE  OF  SHIELD  = 2 
PRESTIGE  ITEMS  = 6 
SCHEMATISM  = 0,1 


ERVIDEL 


0 100 
Output  0 0 10 
0 100 


REGION  = GU  ADI  AN  A 
IMPORTED  ITEMS  = 3 
CHRONOLOGY  = LB  A III 


J.A.  BARCELO,  1995  Back-propagation  algorithms  to  compute  similarity 
relationships  among  archaeological  artifacts.  In  Computer  Applications  in  Archaeology. 


Edited  By  r J.  Wilcock  y K.  Lockyear.  Oxford:  British  Archaeological  Reports. 


DIAZ,D.,  CASTRO, D.,  2001,  “Pattern  Recognition  applied  to  Rock  Art”. 

In  Archaeological  Informatics:  Pushing  the  Envelope.  Edited  by  Goran  Burenhult. 
Oxford:  ArchaeoPress  (BAR  Int.  Series  S1016).,  pp.  463-468. 


LOHSE,  E.S.,  SCHOU,C.,  SCHLADER,R.,  SAMMONS, D,  2004,  “Automated  Classification  of  Stone 
projectile  Points  in  a Neural  Network”.  In  Enter  the  Past.  The  e-way  into  the  four  dimensions  of  culture 
heritage.  Edited  by  Magistrat  der  Stadt  Wien-Referat  Kulturelles  Erbe-Stadtarhchaologie  Wien.  Oxford, 
ArcheoPress  (B  AR  Int.  Series,  S1227),  pp.  431-437). 


I ■ Archaeological  Autoclassification  System 


File  Help  Analysis 
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SIGGI-AACS 


Start  | Load  / Edit  | Threshold]  Image  Smoothing  Results 


Primary  Classification: 

The  Best  Match  for  the  item  is: 
SimpleLanceolate 


Previous  T ab 


Next  T ab 


BP  Pattern  Matching  Demo  - NeuroBox  Neural  Network  Library  Q © 


Network  About 
Input 
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The  general  idea  is 
that  a neural 
network  can  be 
fed  with  visual 
inputs  (images), 
and  it  will  output 
a shape-based 
recognition  of  the 
visual  features 
present  in  that 
input. 


Kashyap,  H.K.,  Bansilal,P.,  Koushik  A.P.,  2003,  Hybrid  Neural  Network  Architecture  for  Age  Identification  of  Ancient  Kannada 
Scripts.  Proceedings  of  the  2003  IEEE  International  Symposium  on  Circuits  and  Systems  ( ISCAS  2003),  Vol.  3,  Page(s):  423-426,  May  25- 
28,  2003. 
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Figure  6.  Set  of  Characters  ranging  join  3rd 
century  BC  to  lS:h  Century  AD. 


An  Artificial  Neural 
Network  for  studying 
ancient  Indian  documents 
written  in  Kannada,  a 
language  of  southern 
India,  which  is  as  old  as 
5th  century  AD. 

(Kashyap  et  al.,  2003). 


Input  Image 


From  Classification  to  Typology 


The  idea  of  Typology.  A set  of  ordered 
Prototypes 


Prototype.  It  is  an  individual  instance  of  some 
entity  serving  as  a typical  example,  for  other 
entities  of  the  same  category.  When  the 
regularities  extracted  for  a given 
archaeological  data  share  a common  set  of 
attributes,  this  set  can  be  said  to  define  a 
prototype. 


From  Classification  to  Typology 


Can  we  define  a “class”  or  “type”  from  a cluster  of  similar  objects? 


In  clustering,  a set  of  explanations  will  be  modeled 
by  first  describing  a set  of  prototypes,  then 
describing  the  objects  using  these  prototypical 
descriptions.  Each  description  gives  the 
probabilities  of  the  observable  features,  assuming 
that  what  has  been  perceived  belongs  to  a group 
composed  of  similar  looking  percepts.  The 
prototype  descriptions  are  chosen  so  that  the 
information  required  to  describe  objects  in  the 
class  is  greatly  reduced  because  they  are  “close”  to 
the  prototype. 


Learning 

We  should  learn  CAUSAL  ASSOCIATIONS. 

B depends  on  A if  B is  consequence  of  A. 

This  relationship  between  EFFECT  and  CAUSE  is  not 
necessary  a formal  relationship.  At  least  not  always.  It 
should  exist  some  type  of  association,  but  we  do  not  need 
that  that  association  follows  formal  criteria.  Obviously,  the 
more  formal  the  dependence  between  EFFECT  and 
CAUSE,  the  greater  validity  we  will  grant  to  the  problem 
solution. 


OBSERVATION  OF  INDIVIDUAL  INSTANCES 


+ 

f>RIOR  KNOWLEDGE: 

constraints  that  will  ensure  that  the  predictions  drawn 
in  automated  archaeologist  will  tend  to  be  plausible  and  relevant  to 

| | system’s  goals  ■ ■ 


testing 

Cause 

Effect 

Feedback  cycle 

Inference  of  a general  model 
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nun  ML 
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Arqueologia  = CIENCIA  EXPERIMENTAL 


FUNCTIONAL  ANALYSIS 


The  function  of  a certain  item  is  - or  should  be  - what  it  is  best  able  to 
do  (or  be)  given  its  physical  constitution  and  its  context. 

D.  Dennet  Design  Stance. 
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FUNCTIONAL  ANALYSIS 


An  object’s  use  can  be  defined  as  the  exertion  of  control  over  a freely 
manipulable  external  object  with  the  specific  intention  of: 

(1)  altering  the  physical  properties  of  another  object, 
substance,  surface  or  medium  via  a dynamic  mechanical 
interaction,  or 

(2)  mediating  the  flow  of  information  between  the  tool  user 

and  the  environment  or  other  organisms  in  the  environment. 
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To  infer  “functionality”  from  perceived 
properties  we  need  a combination  of  three 

kinds  of  information: 


• Knowledge  about  how  the  designers  intended  to  design  the 
artifact  to  have  the  function 

• Knowledge  about  how  the  makers  determined  the  physical 
structure  of  that  artifact  on  the  basis  of  their  technological 
abilities 

• Knowledge  about  how  the  artifact  was  determined  by  its 
physical  structure  to  perform  that  function 
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FUNCTIONAL  ANALYSIS. 

DIRECTLY  INTERACTING  WITH  REAL  OBJECTS 

The  constraints  that  are  relevant  in  function  of  the  archaeological  entity 

fall  into  different  categories,  which  would  include  the  following  (St. 

Amant  2002,  Bicici  and  St.  Amant  2003): 

• Spatial  constraints  describe  the  spatial  relationships  associated  with  a 
tool  and  its  use  in  an  environment.  For  example,  to  use  a hammer  one 
needs  enough  room  to  swing  it. 

• Physical  constraints  describe  physical  relationships  in  the  use  of  the  tool, 
such  as  weight  or  size. 

• Dynamic  constraints  describe  movement-  or  force-related  properties  of 
tool  use.  For  example,  one  needs  to  swing  a hammer  with  appropriate 
speed  in  its  use. 
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The  PEDRA  System 


OUTPUT: 


transversal 


longitudinal 


Original  Image 
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Inputs 
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MEASURES: 

Shape 

Luminosity 


“Texelization” 


DESCRIBING  TEXTURE 


Central  Tendency  texture  values 

In  each  picture  we  calculate  the 
mean  and  standard  deviation 
values  of  texel  features 


transversal 


ANGLE 


longitudinal 
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Mean  = 48.9 
SD  = 30.38 


Mean  = 29 
SD  = 27.01 


NEURAL  NETWORKS 
The  PEDRA  Project 


PEDRA.  ( Stone  in  Catalan) 
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PATTERN  RECOGNITION 


The  network  has  correctly  learn  to  distinguish  use  kynematics 


Edge  orientation 


Edge  orientation 


Hictden  layer 


Original  Image 


Segmentation 


OUTPUT: 


Outputs 

BONE 

SHELL 

FRESH  WOOD 
DRY  WOOD 
FRESH  HIDE 
DRY  HIDE 
MEAT 

\/ 

MEASURES: 

“Texelization” 

Output  / Desired 

BONE 

BUTCHERY 

DRY  HIDE 

HRY  WOOD 

FRESH  HIDE 

FRESH  WOOD 

SHELL 

BONE 

55 

0 

0 

6 

1 

3 

10 

BUTCHERY 

8 

43 

4 

2 

5 

1 

0 

DRY  HIDE 

13 

4 

46 

6 

0 

6 

3 

DRY  WOOD 

13 

1 

12 

43 

1 

13 

4 

FRESH  HIDE 

3 

17 

7 

3 

18 

0 

0 

FRESH  WOOD 

10 

1 

5 

6 

0 

28 

8 

SHELL 

20 

1 

6 

5 

0 

11 

44 

Performance 

BONE 

BUTCHERY 

HIDE 

WOOD 

FRESH  HIDE 

FRESH  WOOD 

SHELL 

MSE 

0,1433 

0,0626 

0,0983 

0,0962 

0,0405 

0,0852 

0,0892 

NMSE 

0,7731 

0,5364 

0,7271 

0,7844 

0,8473 

0,7795 

0,7449 

MAE 

0,237 

0,1455 

0,1966 

0,20543 

0,1137 

0,1910 

0,1950 

Min  Abs  Error 

1, IE-05 

8 , 2E - 05 

0,0006 

0,0006 

0,0002 

2 , 7E - 05 

0,0002 

Max  Abs  Error 

1,0388 

0,9980 

0,9961 

0,9779 

0,9352 

1,0341 

0,9818 

r 

0,512 

0,7116 

0,5397 

0,477 

0,4998 

0,4709 

0,5119 

Percent  Correct 

45,08 

1 

64,1791 

57,5 

60,5633 

72 

45,1612 

63,7681 

BEYOND  VISUAL  DATA 


Download  from 

Dreamstime.com 

This  watermarked  comp  image  is  for  previewing  purposes  only. 
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John  Takai  | Dreamstime.com 


127162 


Intensity 


BEYOND  VISUAL  DATA: 


Compositional  Data 


Li 


Na 


1 1 

u 

m, 

1 3- 

a 

. a 

k 

Na 


La 


CaO 


Sb 


Fb 


iiiSfo 


r 


k 2Sn 
PhKIiOj 


FI 

41 

* 

1 

Y 

Sh 

i 7 

. ,i.  .i.  . 

VUtO 

n 


ICO 


200 


300 


-rjjjl 


400 


>l;blri\ 


u 


Ofc 


/ 


r 


* *0 


500 


-i^i L ^ ^-|L,L  -v  » ,J 1. 


600  900  1000  ItTO  1200  1300  1400  1500  1000  1700  1600  1900  2000 

m iz 

127  7163 


100 


£00 


300 


400  Mass  (rntej  | 


127/ 


Mechanical 

properties 

of 

archaeological 

materials 


164 


Prfldicpl  rrcililF-m 


CONSTRAIN  1$ 


Task  <.'i>rit  trjijili 
{ Lewd  oJ  Acceptability} 

Task  hfatambn 

- f recauun 

^ ■ Fa  ret  — 

- NftlU'f  of  Ac4HHi 

BnfleAcacy 

Quantity 

Time  Available 

failure  Cunsequerura  (Rk1c) 

Laedlicnil  Constrainls, 

(WKHiaiyw  siirta^k  l«*liCUi&  Or  flWirtminenb} 

Material  Constraints 

Available  Md  ri-sls,  5izct,  and  CutitE 
[irLlivr  Performance 
Sdil  ivu  Weai/Failure  Kales 
raseffTWork^ 

TednrttUigical  Gmitrainffs 
Available  Techiiahsgy 
IWucHonCo&ls 

R'Opdir/E'Uf^idrpfflinfj/ElepLannieni  CohI  . 
bkilJ  Required. 


Sac  :occ  uiiuni  i ; CanettilolS 

Mobility 

Ittiispon  Capacity 
Available  Libor 
Siic>ragr/uw#n|cw^  capability 


Ideologic' I Lonatrdinb 


Design  Cvraiilrralinna 
-Size  arid  Weight 

- Edge  A i'-K.li!  and  Farm 

- Prehension  and  Hefting 

- Length  of  Use  (Lbe-lile)  — ► 

« Dc'^rt1:  ufJjptjalrd  Use 

-Reliability 

frubuvlrietib  and  ‘rjvefdesi^n") 
i liaie  ijl  ri,p,in1'ri  -,hai  rv  ii.ind;  iftivirtmeMbi 
-MnltifuflCtKiridJity  (vccauhltty) 


Trad  ltd  ioji/Reduclian 
- frfeeh'rpcrtmg  Strategies 

eg  - Expedient  Slock  Cose 
t Bifore 

- [.cing'  [.j  iifdi 

flak*/ Blade- Ttoli 

- scuviMifl.r.g/  Etarycling 
» Bipolar 

- GmuutriflAe 

• Etaharpming: 

1 hnrd  hummer  (notching, 
ttfilinunus  ruiuudung, 
burinatui^.  etc.) 

-tilLlei 


- pruauni 
• glUVlLfVg 


A schematic  represmtatioB  of  ihe  design  arri  jModucriwi  piwcsj  for  prwIkaJ  lithic  technology.  Tk  \a  a se- 

quential design  process  tbai  rciulLi  In  [be  pccxtucLiu^ii' s tool.  Once  it  is  completed  and  the  lout  is  cried  out,  b may  lead  io 
evalujdtwi  uf  Ihe  aL-ueplabilily  «F  results  and  Tepswssineftl  of  various  factor,  Leading.  EO  a new  inund  n\  tnnl  pncschidinn. 


FUNCTIONAL  ANALYSIS. 

DIRECTLY  INTERACTING  WITH  REAL  OBJECTS 


MEASURING  MECHANICAL  PROPERTIES 


igure  I . (a}  The  model  jar  manufactured  for  the  purposes  of  developing  and  testing  the  FEA  method,  (b)  Testing  configuration  of  the  test  jar. 


3 - point 
Flexure  test 
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MEASURING  PHYSICAL  PROPERTIES 

Physical  properties  - are  those  whose  particular  values  can  be  determined 
without  changing  the  identity  of  the  substance:  DENSITY,  MOISTURE 
CONTENT,  PERMEABILITY,  SHRINKAGE 

Friction  properties  - include  the  coefficients  of  static,  kinetic,  and  rolling  friction, 

Thermal  properties  - thermal  conductivity,  thermal  diffusivity,  thermal  expansion 
coefficient,  thermal  shock  resistance,  specific  heat,  melting  point,  creep  resistance. 
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MEASURING  MECHANICAL  PROPERTIES 

Mechanical  properties  - ThE  value  may  vary  as  a result  of  the  physical 
properties  inherent  to  each  material,  describing  how  it  will  react  to  physical 
forces.  The  main  characteristics  are  ELASTIC,  STRENGTH  and  VIBRATION. 

ELASTIC  PROPERTIES : 

modulus  of  elasticity,  is  the  ratio  of  linear  stress  to  linear  strain. 

Poisson's  ratio  is  the  ratio  of  lateral  strain  to  axial  strain. 

Yield  strength  refers  to  the  point  on  the  stress-strain  curve  beyond  which  the  solid  starts 
to  deform  plastically  and  cannot  be  reversed  upon  removal  of  the  Loading 


127/  168 


MEASURING  MECHANICAL  PROPERTIES 

STRENGTH  PROPERTIES:  The  material’s  mechanical  strength  properties  refer 
to  the  ability  to  withstand  an  applied  stress  without  failure,  by  measuring  the 
extent  of  a material's  elastic  range,  or  elastic  and  plastic  ranges  together.  Loading, 
which  refers  to  the  applied  force  to  an  object,  can  be  by: 

Tension  -It  can  be  quantified  as  uWmate  tensile  strength,  which  is  the  maximum 
amount  of  tensile  stress  a material  can  withstand  while  being  stretched  or  pulled 
before  failure.  Ductility  measures  how  much  a material  deforms  under  tensile 
load  before  breaking.  It  can  be  measured  in  percentage  of  elongation  of  a tensile 
sample  after  breaking.  On  the  contrary,  brittleness  is  the  ability  of  a material  to 
fracture  with  very  little  or  no  previous  detectable  deformation. 
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MEASURING  MECHANICAL  PROPERTIES 

Compression  - involves  pressing  the  material  together.  In  fact,  it  is  the  opposite 
of  tensile  loading.  Compressive  strength:  is  the  maximum  amount  of  compressive 
stress  a material  can  withstand  while  being  compressed  before  failure. 

Hardness  :A  measure  for  material  hardness  can  also  be  the  degree  of  abrasion, 
which  is  the  resistance  to  grinding  force. 

Bending  - involves  applying  a load  that  causes  a material  to  curve,  resulting  in 
compressing  the  material  on  one  side  and  stretching  it  on  the  other.  It  can  be 
quantified  as  bending  strength  and  flexural  strength. 
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FUNCTIONAL  ANALYSIS. 

DIRECTLY  INTERACTING  WITH  REAL  OBJECTS 
MEASURING  MECHANICAL  PROPERTIES 

Shear  - involves  applying  a load  parallel  to  a plane,  causing  the  material  on  one 
of  the  sides  of  the  plane  to  want  to  slide  across  the  material  on  the  other  side.  It 
can  be  quantified  as  Shear  strength,  which  is  the  maximum  amount  of  shear  stress 
a material  can  withstand  before  failure.  Shear  strain:  change  in  the  angle  between 
two  perpendicular  lines  in  a plane.  Shear  modulus  (or  modulus  of  rigidity,  ratio 
of  shear  stress  to  shear  strain,  measures  the  stiffness  of  materials  indicating  the 
resistance  to  deflection  of  a member  caused  by  shear  stresses. 

Torsion  - Torsion  strength  indicates  the  applied  force  which  causes  twisting  in  a 
material. 

Fatigue  - Fatigue  IJmit  refers  to  the  maximum  stress  a material  can  withstand 
under  cyclic  loading.  This  resistance  to  failure  under  particular  combinations  of 
repeated  loading  conditions  is  measured. 
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FUNCTIONAL  ANALYSIS. 

DIRECTLY  INTERACTING  WITH  REAL  OBJECTS 


Imagine  the  answer  of  a Museum  director  when  we  ask  her  to  break  a prehistoric  object 
in  order  to  measure  its  structural  and  mechanical  properties.  Given  that  prehistoric  and 
ancient  objects  tools  not  always  can  be  used  in  the  present  nor  “touched”  to  preserve  its 
integrality,  we  are  limited  to  the  possibily  of  manipulating  a virtual  surrogate  of  the 
Object.  127/172 
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Reverse  Engineering  is  the  process  of  extracting  missing  knowledge  from 
anything  man-made,  by  going  backwards  through  its  development  cycle  and 
analyzing  its  structure,  function  and  operation.  It  consists  of  a series  of  iterative 
steps,  each  addressing  different  questions  regarding,  in  this  case,  an  overall 
artefact.  These  steps  may  be  repeated  as  often  as  needed  until  all  steps  are 
sufficiently  satisfied. 
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FUNCTIONAL  ANALYSIS. 
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FUNCTIONAL  ANALYSIS. 

SIMULATING  THE  INTERACTION  WITH  SURROGATES  OF 

REAL  OBJECTS 


We  can  imagine  different  kinds  of  simulated 

interaction: 

1 ) static,  which  calculates  displacements, 
reaction  forces,  strains,  stresses,  and 
factor  of  safety  distribution; 

2)  frequency,  calculates  stresses  caused  by 
resonance; 

3)  fatigue,  calculates  the  total  lifetime, 
damage,  and  load  factors  due  to  cyclic 
loading; 

4)  non-linear,  calculates  displacements, 
reaction  forces,  strains,  and  stresses  at 
incrementally  varying  levels  of  loads  and 
restraints; 

5)  dynamic,  calculates  the  model's  response 
due  to  loads  that  are  applied  suddenly  or 
change  with  time  or  frequency 


FUNCTIONAL  ANALYSIS. 

SIMULATING  THE  INTERACTION  WITH  SURROGATES  OF 

REAL  OBJECTS 

J.A.  BARCELO  & V.  MOITINHO  2011-2013:  Prehistoric  Bow  and  Arrow 
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FUNCTIONAL  ANALYSIS. 

SIMULATING  THE  INTERACTION  WITH  SURROGATES  OF 

REAL  OBJECTS 

La  Draga  (Banyoles).  An  early  neolithic  site 


Bow  (D/02 -]] 89 -11), 

3D  Scans  alignment. 

Total  number  of  scans:  179,  in  3 
parts. 

FOV:  90  mm,  stereo. 

Resolution  50  jum. 


Bow  (D/02 -JJ 89-11). 

Merged  polygonal  mesh  - 3D 
digital  surface  model. 


Polygonal  mesh  with  poly-faces  , 
detail. 
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FUNCTIONAL  ANALYSIS. 

SIMULATING  THE  INTERACTION  WITH  SURROGATES  OF 

REAL  OBJECTS 

.A.  BARCELO  & V.  MOITINHO  2011-2013:  Prehistoric  Bow  and  Arrow 


[M]{a}  + [C]M  + [K]{x}  = {R 
Where: 

M = Mass  matrix  a = acceleration  vector 


C = damping  matrix  v = velocity  vector 

K = Stiffness  matrix  x = displacement  vector 

F=  external  force  vector 


Pressure,  displacement  and  deformation 
analysis. 
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